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1. Introduction 

There has been a steady progression of lattice results from a quenched sea to two-flavour 
and more recently 2+1 flavour sea in an attempt to provide a more complete and quantitative 
description of hadronic phenomena. (By 2 + 1 flavours we mean here 2 mass degenerate up-down, 
mi, quarks and one strange, m^, quark.) In this talk we shall consider a 2 + 1 flavour clover action 
and discuss some ways of approaching in the mi-nis plane the physical point {m'^,m*), where the 
natural starting point for these paths is an SUf{3) flavour symmetric point m/ = mv = mfym. For 
clover (i.e. non-chiral) fermions a problem arises because the measured pseudoscalar masses are 
proportional to the renormalised quark masses. The singlet, S, and nonsinglet, NS quark masses 
renormalise differently which means that the relation to the bare quark masses and hence k, which 
is the adjustable simulation parameter, is more complicated Choosing the path such that the 
singlet quark mass is kept fixed provides an elegant solution to this problem. This procedure 
also has the advantage that it enables a wide range of quark masses to be probed (including the 
mass of the strange quark) and is thus particularly useful for strange quark physics and SUf{3)- 
chiral perturbation theory, as the kaon mass is never larger than its physical value. (Note that 
this path choice although favourable for clover-type fermions is not restricted to them.) Indeed 
this includes the case with a heavy light quark mass and hght strange quark mass. As a 'proof 
of concept' preliminary results given here show the correct splittings for the baryon (octet and) 
decuplet spectrum. 

The particular clover action used here has a single iterated mild stout smearing for the hopping 
terms together with thin links for the clover term (this is in an attempt to ensure that the fermion 
matrix does not become too extended). Together with the (tree level) Symanzik improved gluon 
action this constitutes the Stout Link Non-perturbative Clover or SLiNC action for which the clover 
coefficient, Csw, has recently been non-perturbatively (NP) determined, [^, using the Schrodinger 
Functional, or SF, formalism. Further details about the action may be found in this reference. 
Simulations have been performed using HMC with mass preconditioning for 2 mass-degenerate 
flavours and the rational approximation for the 1-flavour. Two programmes were used, a Fortran 
programme, [^, and also the Chroma programme, [Q]. Quark mass degenerate runs (denoted by the 
subscript 'sym') on 16'' x 32 lattices followed by 24^ x 48 lattices have located a suitable j8-range, 
see Fig. |. We see that there is good agreement between SF and pseudoscalar mass determinations 
of the critical value of Kgym, ^sym;c- The simulations reported here have been performed at j8 = 5.50 
which gives a ~ 0.08 fm (taking the scale ro = 0.5 fm and a linear extrapolation of ro/a to the 
chiral limit). We have furthermore checked that there is a distinct gap between the distribution of 
(the modulus of) the lowest eigenvalue and indicating that the simulations are stable on present 
volumes. Further details of the results in this write-up will be given in [^. 

2. Non-degenerate quark masses 

As mentioned before the problem (at least for clover-like fermions, with no chiral symmetry) 
is that singlet and non-singlet quark mass can renormalise differently. 
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Figure 1: Results for the pseudoscalar mass against 1 / (fv,.,,, for j3 = 5.40, 5.50, 5.60. For j3 = 5.50 a linear 
fit (line) using the lightest three masses is made. For comparison a quadratic fit (dashed line) is also made to 
all the masses. The stars represent the SF determination of Ksym;c\ the open circle is the result from the linear 
fit. 
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(2.1) 



with q G {sym,l,s,v} (also including possible different valence quarks to sea quarks) and 

rn = -{2mi +ms) . 

(We are assuming here that the quadratic improvement terms, [^], are small, [^.) The bare quark 



(2.2) 



mass in eq. (2. 1) is defined by 



1 f 1 



ama 



(2.3) 



Now if nips {I, s) is the measured pseudoscalar mass (with quarks q\, ^2) then we expect that 



In particular 



(anij^)^ = (am^^)^ oc 2(am^ + azam) 



(2.4) 



(2.5) 



which is <^ aniq unless q = I = s = sym (i.e. mq = m). That az is non-zero may be easily seen 
by considering partially quenched results (i.e. results for the pseudoscalar mass where the valence 
quark masses may differ from the sea quark masses). In Fig. |2| we show some partially quenched 
results and compare the extrapolation with the 5?7/r(3)-symmetric results. The lines clearly have 
different gradients. OCz can be estimated as {aruj^Y vanishes at say {pq for 'partially quenched'), 
giving 
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Figure 2: The left plot shows the Slip (3) symmetric pseudoscalar masses versus anig (in red) together with 
the partially quenched results (in black) using the mass degenerate sea quark fc,,™ — 0.12090. (This is the 
kappa value at the crossing point.) The right plot shows equivalent domain wall results from JtI table V. 



(setting I = s = sym in m). This gives here az ~ 1-2 (but the determination is quite sensitive to 
small changes in Ksym-^c and K^^). This is to be compared with domain-wall fermions, also shown 
in Fig. ^ where the results line up. 

We wish to approach the physical point along some path in the mf-mf plane (m^ is considered 
as it is related to the measurable pseudoscalar mass) from an SUf{3) symmetric point (mf = mf = 
nhym)- This is depicted in Fig. ^ Two possibilities are mf = const. = mf * (i.e. strange quark mass 




Figure 3: The left sketch shows the mf-mf plane. The physical quark masses are denoted by (mf*,mf*). 
The dashed diagonal Une is the 5t/F(3)-symmetric line. Two possible paths from this line are shown, mf = 
mf * and m^ — const. . The right sketch shows the equivalent results in the mi-nis plane. 

being held constant) or m^ = const, (i.e. singlet quark mass being held constant). Note that the 
region covered is mf > 0, mf > 0. For simulations we need to translate this to unrenormalised 
quantities also shown in Fig. |3[ Note that the physical domain is now m/ > — (|az/(l + \oiz))ms, 
rus > — (|az/(l + lOJz))"^/- While mf = const. = mf * translated to bare quark masses now depends 
on the difficult-to-determine az, the singlet quark mass rn'^ oc m or m^ = {2m"^ + m*) — 2m/ is 
independent of the value of az- This motivates the choice of this path. As m/ \ m^* then m^ ^ m* 
i.e. the m^-m/ splitting or mk increases to its physical value. Other potential advantages include: 
the singlet quark mass is correct from the very beginning; flavour singlet quantities are flat at the 
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symmetric point - allowing simpler extrapolations; numerically the HMC cost change should be 
moderate along this path. (These points will be further discussed in |^.) 

Of course practically we must now determine the initial Ksym to complete the relation between 
Ks and Ki, 

K-. = ^^. (2.7) 

„(0) K, 



For our path choice, K^^m can be implicitly found by using eqs. {lA), (2S) together with a singlet 



scale X (so X* = X^yli here) to relate the known physical point to the initial symmetric point, 



1 _ \{2ml + ml) 



_R(0) a< ,^ 
oc m = msym °^ — . (2.8) 



cx " 

For X we have several choices. For example 

• The centre of mass^ of the octet baryons: = | (m^ + m| + mi ) (stable under strong interac- 

(0)^2/^^^(0)^2 



tion) giving i(2m| + m|)/i(m^ + m| + m|)|* = 0.169/1.34 = 1/7.93 = (amD7(am 



N 



Centre of mass of the decuplet baryons: X^ = \ (2m^ (which decay under strong inter- 
actions) giving i(2m| + m2)/l(2m| + m2j)|* =0.169/1.94 = 1/11.5 = {am'^^^Y /{am^Y, 

Agluonic quantity X2 = l/r^ so r2(2m| + m|)|* = 0.169/0.395^ = 1.083 = {rf^ /af{amfY 
(ro =0.5fm). 



Thus where the numerically determined 5'?7/r(3)-symmetric line crosses with the line {aX^ 



\2 



sym ) 



cxiamnf' gives our initial point. The results are shown in Fig. ^ We see that they are all consistent 
around the lightest pseudoscalar mass, namely for = 0.12090 which we shall take as our start- 



ing value. From eq. ( |2.7[ ) we now have a relation between and Ky. After some experimentation 
we chose the Ki, values given in Table |I| Note that is possible to choose K;, Kg values (here 



0.12083 0.12104 mi > m, 

0.12090 0.12090 m, =m, 

0.12095 0.12080 m/ < m, 

0.12100 0.12070 mi<m, 

0.12104 0.12062 mi<ms 

Table 1: Present (fC/, k^) values (simulated on 24^ x 48 lattices). 

(0.12083,0.12104)) such that m/ > m,. In this strange world we would expect to see an inversion 
of the particle spectrum, with, for example, the nucleon being the heaviest octet particle. 

3. Hadron spectrum 

As an example we now give some results in Fig. |5| for the baryon decuplet: Q., S*, £*, A. 
These are at present low statistics results and are meant to illustrate the 'proof of concept' only. 
Noteworthy is that the correct ordering of the particle spectrum has been achieved. 
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Figure 4: {arriN)^ against {amnY (upper left picture) together with the hne {amt^Y = 7 .93{amji)^; (arri/^)^ 
against iam^)^ (upper right picture) together with the line {am\)^ = \ \.5{amTc)^; ro/a against [anijif' (lower 
picture) together with the line (ro / a)^ = 1 .083 / {am;^)^. All data comes from the symmetric points. 



4. Conclusions 



For the NP O (a) -improved 2+1 flavour clover action discussed here we have first found that 
there is consistency for the critical k;„,„, between the SF approach and the determination from the 
vanishing of the hadron pseudoscalar mass. The quark mass renormalisation suggests that the 
simplest way of approaching the physical point is to hold the singlet quark mass fixed. Exploratory 
results for the hadron mass spectrum give splittings in the correct order (including inversion when 
mi > nis i.e. we can simulate a strange world where, for example, the nucleon can decay). Further 
results will be published elsewhere, [§]. 
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Figure 5: The decuplet baryons. The upper left picture shows the A, upper right E*, lower left E* and lower 
right the £2. The common SUp (3) symmetric quark mass values are shown in black. The experimental values 
are shown with stars. 
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